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ABSTRACT 

A century ago, in 1911 and 1913, Plummer and then Reynolds introduced their models to 
describe the radial distribution of stars in 'nebulae'. This article reviews the progress since 
then, providing both an historical perspective and a contemporary review of the stellar structure 
of bulges, discs and elliptical galaxies. The quantification of galaxy nuclei, such as central mass 
deficits and excess nuclear light, plus the structure of dark matter halos and cD galaxy envelopes, 
are discussed. Issues pertaining to spiral galaxies including dust, bulge-to-disc ratios, bulgeless 
galaxies, bars and the identification of pseudobulges are also reviewed. An array of modern 
scaling relations involving sizes, luminosities, surface brightnesses and stellar concentrations are 
presented, many of which are shown to be curved. These 'redshift zero' relations not only quantify 
the behavior and nature of galaxies in the Universe today, but are the modern benchmark for 
evolutionary studies of galaxies, whether based on observations, A-body-simulations or semi- 
analytical modelling. For example, it is shown that some of the recently discovered compact 
elliptical galaxies at 1.5 < z < 2.5 may be the bulges of modern disc galaxies. 

Subject headings: galaxy structure — galaxy scaling relations — galaxy elliptical — galaxy spiral — 
galaxy compact — galaxy cD, halos — galaxy nuclei — galaxy central mass deficits — galaxy excess 
nuclear light — galaxy bulge-disc ratios — galaxy bars — galaxy pseudobulges — galaxy bulgeless — 
galaxy dust — Sersic model — core-Sersic model — Einasto model — dark matter halos. 



1. Introduction 

For the last century astronomers have been 
modelling the structure of 'nebulae', and here we 
focus on those external to the Milky Way. A 
key activity performed by many astronomers, past 
and present, is the catergorisation of these galax- 
ies (Sandage 2005) and the quantification of their 
physical properties. How big are they? How 
bright are they? What characteristics distinguish 
or unite apparent subpopulations? Answers to 
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such questions, and the establishment of "scaling 
relations" between two or more galactic properties 
provides valuable insight into the physical mecha- 
nisms that have shaped galaxies. 

Understanding how galaxies form, increasingly 
through the use of simulations and semi-analytic 
modelling (e.g. Cole 1991; White & Frenk 1991; 
Kauffmann et al. 1993, 2003; Avila-Reese et al. 
1998; Cole et al. 2000; de Lucia et al. 2006; Bower 
et al. 2006; Kauffmann et al. 2004; Di Matteo et 
al. 2005; Croton et al. 2006; Naab et al. 2006; 
Nipoti et al. 2006; Covington et al. 2011; Guo et 
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al. 2011, and references therein), requires an ac- 
curate knowledge of galaxy properties and scaling 
laws, as elucidated by Driver (2011). Not surpris- 
ing, our knowledge of galaxies is best in the nearby 
Universe — out to distances typically measured in 
megaparsecs rather than by redshift z — where 
galaxy structures can be reasonably well resolved. 
The properties of these galaxies provide the z = 
benchmark used in the calibration of galactic evo- 
lutionary studies — both observed and simulated. 

Popular scaling relations involving global galaxy 
parameters such as size, surface brightness, lumi- 
nosity and concentration are reviewed here. As we 
shall see, many bivariate distributions, which are 
frequently assumed to be linear, are often only ap- 
proximately so over a restricted luminosity range. 
For example, it may come as a surprise for many 
too learn that the useful Kormendy (1977b) re- 
lation is only the tangent to the bright arm of 
a continuous but curved effective radius- (surface 
brightness) relation which unifies dwarf and gi- 
ant elliptical galaxies (section 13.2.41) . Similarly, 
the Faber- Jackson (1976) relation with a slope of 
4 represents the average slope over a restricted 
luminosity range to what is a curved or broken 
luminosity- (velocity dispersion) distribution, in 
which the slope is 2 rather than 4 at lower lumi- 
nosities (section [333]). Knowing these trends, the 
bulk of which cannot be established when assum- 
ing structural homology, i.e. using de Vaucouleurs' 
(1948) R 1 / 4 model, is vital if one is to measure, 
model and make sense of galaxies. 

This article has been structured into four main 
sections. Section 1 provides this general overview 
plus a further review and introduction to galax- 
ies on the Hubble- Jeans sequencfl Included are 
diagrams showing the location of dynamically hot 
stellar systems in the mass-size and mass-density 
plane, revealing that some high- z compact galax- 
ies have properties equivalent to the bulges of lo- 
cal disc galaxies. Section [5] provides an histori- 
cal account of how the radial distribution of stars 
in elliptical galaxies have been modelled, and the 
iterative steps leading to the development of the 



1 This review does not encompass dwarf spheroidals, or any, 
galaxies fainter than Mb ~ —14 mag. These galaxies can 
not be observed (to date) at cosmologically interesting dis- 
tances, and their increased scatter in the colour-magnitude 
relation may indicate a range of galaxy types (e.g. Penny 
& Conselice 2008, and references therein). 



modern core-Sersic model ( section |2.2[) . Subsec- 
tions cover the Sersic model (section |2.1|) . its re- 
lation and applicability to dark matter halos (sec- 
tion 12.1.11) . partially-depleted galaxy cores (sec- 
tion 12.2. IT) , excess nuclear light (section |2"3|) and 
excess light at large radii in the form of halos 
or envelopes around giant elliptical galaxies (sec- 
tion [231) • Section [3] presents and derives a number 
of elliptical galaxy scaling relations pertaining to 
the main body of the galaxy. From just two linear 
relations which unite the faint and bright ellipti- 
cal galaxy population (section I3.1[) , a number of 
curved relations are derived (section I3.2[) . Sev- 
eral broken relations, at Mb ~ —20.5 mag, are 
additionally presented in section 13.31 For those 
interested in a broader or different overview of 
elliptical galaxies, some recent good reviews in- 
clude Renzini (2006), Cecil & Rose (2007), Ciotti 
(2009) and Lisker (2009). Finally, the latter third 
of this paper is tied up in section 0] which con- 
tains a discussion of the light profiles of disc galax- 
ies and their bulge-disc decomposition (14.11) . Also 
included are subsections pertaining to dust (sec- 
tion 14. 2j) . the difficulties with identifying pseu- 
dobulges (section |4~3|) . potential bulgeless galax- 
ies (section dm and methods to model bars (sec- 
tion 14. 5|) . Throughout the article references to 
often overlooked discovery or pioneer papers are 
provided. 

1.1. Early Beginnings 

Looking out into the Milky Way arced across 
our night sky, the notion that we are residents 
within a pancake-shaped galaxy seems reasonable 
to embrace. Indeed, back in 1750 Thomas Wright 
also conjectured that we reside within a flat layer 
of stars which is gravitationally bound and rotat- 
ing about some centre of mass. However, analo- 
gous to the rings of Saturn, he entertained the idea 
that the Milky Way is comprised of a large annu- 
lus of stars rotating about a distant centre, or that 
we are located in a large thin spherical shell rotat- 
ing about some divine centre (one of the galactic 
poles). While he had the global geometry wrong, 
he was perhaps the first to speculate that faint, ex- 
tended nebulae in the heavens are distant galaxies 
with their own (divine) centers. 

As elucidated by Hoskin (1970), it was Im- 
manuel Kant (1755), aware of the elliptically- 
shaped nebulae observed by Maupertuis, and 
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working from an incomplete summary of Wright 
(1750) that had been published in a Hamburg 
Journal, who effectively introduced the modern 
concept of disc-like galactic distributions of stars 
— mistakenly crediting Wright for the idea. 

Using his 1.83 m "Leviathan of Parsonstown" 
metal reflector telescope in Ireland, Lord William 
Henry Parsons, the 3rd Earl of Rosse, discovered 
226 New General Catalogued (NGC: Dreyer 1888) 
and 7 Index Catalogue (IC: Dreyer 1895, 1908) ob- 
jects (Parsons 1878). Important among these was 
his detection of spiral structure in many galaxies, 
such as M51 which affectionately became known 
as the whirlpool galaxy. 

Further divisions into disc (spiral) and elliptical 
galaxy types followed (e.g. Wolf 1908; Knox Shaw 
1915; Curtis 1918; Reynolds 1920; Hubble 1926)0 
and Shapley & Swope (1924) and Shapley (1928) 
successfully identified our own Galaxy's (gravita- 
tional) center towards the constellation Sagittarius 
(see also Seares 1928). 

With the discovery that our Universe contains 
Doppler shifted 'nebulae' that are expanding away 
from us (de Sitter 1917; Slipher 1917; see also 
Friedmann 1922, Lundmark 1924 and the reviews 
by Kragh & Smith 2003 and Shaviv 2011), in ac- 
cord with a redshift-distance relation (Lemaitre 
1927, Robertson 1928, Humasson 1929, Hubble 
1929)H — i.e. awareness that some of the "neb- 
uale" are external objects to our galaxy — came 
increased efforts to catergorise and organise these 
different types of "galaxy". As noted by Sandage 
(2004, 2005), Sir James Jeans (1928) was the first 
to present the (tuning fork)-shaped diagram that 
encapsulated Hubble's (1926) early-to-late type 
galaxy sequence, a sequence which had been in- 



2 Freye Urtheile, Achtes Jahr (Hamburg, 1751), translated 

by Hastie, op. cit., Appendix B. 
3 The NGC built upon the (Herschel family's) Catalog of 

Nebulae and Clusters of Stars (Herschel 1864). 
4 Reynolds (1927) called Hubble's attention to pre-existing 

and partly-similar galaxy classification schemes which were 

not cited. 

5 It is of interest to note that Hubble (1934, 1936b, 1937) was 
actually cautious to accept that the redshifts corresponded 
to real velocities and thus an expanding Universe as first 
suggested by others. He used the term "apparent velocity" 
to flag his skepticism. In point of fact, Hubble & Tolman 
(1935) wrote that the data is "not yet sufficient to permit 
a decision between recessional or other causes for the red- 
shift" . 



spired in part by Jeans (1919) and later popu- 
larised by Hubble (1936a; see Block et al. 2004). 
Quantifying the physical properties of galaxies 
along this sequence, with increasing accuracy and 
level of detail, has occupied many astronomers 
since. Indeed, this review addresses aspects re- 
lated to the radial concentration of stars in the 
elliptical and disc galaxies which effectively define 
the Hubble- Jeans sequence. Irregular galaxies are 
not discussed here. 

1.2. The modern galaxy 

For reasons that will become apparent, this re- 
view uses the galaxy notation of Alan Sandage 
and Bruno Binggeli, in which dwarf elliptical (dE) 
galaxies are the faint extension of ordinary and 
luminous elliptical (E) galaxies, and the dwarf 
spheroidal (dSph) galaxies — prevalent in our Lo- 
cal Group (Grebel 2001) — are found at mag- 
nitudes fainter than Mb ~ —13 to —14 mag (« 
1O 8 M in stellar mass; see Figure QJa,). Figure QJi 
reveals a second branch of elliptically-shaped ob- 
ject stretching from the bulges of disc galaxies and 
compact elliptical (cE) galaxies to ultra compact 
dwarf (UCD) objects (Hilkcr ct al. 1999; Drinkwa- 
ter et al. 2000; Norris & Kannappan 2011 and ref- 
erences therein). A possible connection is based 
upon the stripping of a disc galaxy's outer disc 
to form a cE galaxy (Nieto 1990; Bekki et al. 
2001b; Graham 2002; Chilingarian et al. 2009) 
and through greater stripping of the bulge to form 
a UCD (Zinnecker et al. 1988; Freeman 1990; 
Bassino et al. 1994; Bekki 2001a). It is thought 
that nucleated dwarf elliptical galaxies may also 
experience this stripping process, giving rise to 
UCDs. 

While the identification of local spiral galax- 
ies is relatively free from debate, the situation is 
not so clear in regard to elliptically-shaped galax- 
ies. The discovery of UCDs, which have sizes 
and fluxes intermediate between those of galax- 
ies and (i) the nuclear star clusters found at the 
centres of galaxies and (ii) globular clusters (GCs: 
e.g. Hasegan et al. 2005; Brodie & Strader 2006), 
led Forbes & Kroupa (2011) to try and provide a 
modern definition for what is a galaxy (see also 
Tollerud et al. 2011). Only a few years ago there 
was something of a divide between GCs and UCDs 
— all of which had sizes less than ~30 pc — and 
galaxies with sizes greater than 120 pc (Gilmore et 
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Fig. 1. — Left panel: The radius containing half of each object's light, R1/2 (as seen in projection on the sky), 
is plotted against each object's stellar mass. Open circles: dwarf elliptical (dE) and ordinary elliptical (E) 
galaxies from Binggeli & Jerjen (1998), Caon et al. (1993), D'Onofrio et al. (1994) and Forbes et al. (2008). 
Filled circles: Bulges of disc galaxies from Graham & Worley (2008). Shaded regions adapted from Binggeli 
et al. (1984, their figure 7), Dabringhausen et al. (2008, their figure 2), Forbes et al. (2008, their figure 7), 
Misgcld & Hilker (2011, their figure 1). The location of the so-called "compact elliptical" (cE) galaxies is 
shown by the rhombus overlapping with small bulges. The location of dense, compact, z = 1.5 galaxies, as 
indicated by Damjanov et al. (2009, their figure 5), is denoted by the dashed boundary overlapping with 
luminous bulges. Right panel: Stellar mass density within the volume containing half each object's light, 
rhoi/2, versus stellar mass. The radius of this volume was taken to equal 4/3 xR 1 / 2 (Ciotti 1991; Wolf et 
al. 2010). 
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al. 2007). However, as we have steadily increased 
our celestial inventory, objects of an intermediate 
nature have been found (e.g. Ma et al. 2007, their 
Table 3), raising the question asked by Forbes & 
Kroupa for which, perhaps not surprisingly, no 
clear answer has yet emerged. While those au- 
thors explored the notion of a division by, among 
other properties, size and luminosity, they did not 
discuss how the density varies. As an addendum 
of sorts to Forbes & Kroupa (2011), the density 
of elliptically-shaped objects is presented here in 
Figure QJ>. This is also done to allow the author 
to wave the following flag. 

Apparent in Figure [T|d, but apparently not 
well recognised within the community, is that 
the bulges of disc galaxies can be much denser 
than elliptical galaxies. If the common idea of 
galaxy growth via the accretion of a disc, per- 
haps from cold-mode accretion streams, around 
a pre-existing spheroid is correct (e.g. Navarro & 
Bcnz 1991; Steinmetz & Navarro 2002; Birnboim 
& Dekel 2003; see also Conselice et al. 2011 and 
Pichon et al. 2011), then one should expect to find 
dense spheroids at high-z with 10 10 -10 11 M Q of 
stellar material, possibly surrounded by a faint 
(exponential) disc which is under development. 
It is noted here that the dense, compact early- 
type galaxies recently found at redshifts of 1.4-2.5 
(Daddi et al. 2005; Trujillo et al. 2006) display 
substantial overlap with the location of present 
day bulges in Figure QJ,, and that the merger sce- 
narios for converting these compact high-z galax- 
ies into today's elliptical galaxies are not with- 
out problems (e.g. Nipoti et al. 2009; Nair et al. 
2011). It is also noted that well-developed discs 
and disc galaxies are rare at the redshifts where 
these compact objects have been observed along- 
side normal-sized elliptical galaxies. Before trying 
to understand galaxy structure at high-redshift, 
and galaxy evolution — themes not detailed in 
this review — it is important to first appreciate 
galaxy structures at z — where observations are 
easier and local benchmark scaling relations have 
been established. 

2. Elliptical Galaxy Light Profiles 

Over the years a number of mathematical func- 
tions have been used to represent the radial distri- 
bution of stellar light in elliptical galaxies, i.e. their 



light profiles. Before getting to de Vaucouleurs' 
i? 1 / 4 model in the following paragraph, it seems 
apt to first quickly mention some early competi- 
tors. Although Plummer's (1911) internal-density 
model was developed for the nebulae which be- 
came know as globular clusters, because of its sim- 
plicity it is still used today by some researchers to 
simulate elliptical galaxies, even though, it should 
be noted, no modern observers use this model to 
describe the radial distribution of light in ellip- 
tical galaxies. Reynold's (1913) surface-density 
model, sometimes referred to as Hubble's (1930) 
model or the Reynold-Hubble model, was used 
to describe the nebula which became known as 
elliptical galaxies. It has an infinite mass and 
is also no longer used by observers today. The 
modified Hubble model (Rood et al. 1972), which 
also has an infinite mass, is also still sometimes 
used by simulators, even though, again, observers 
do not use this model anymore. Oemler's (1976) 
exponentially-truncated Hubble model, known as 
the Oemler-Hubble model, is also not used to rep- 
resent the observed stellar distribution in elliptical 
galaxies because it too, like its predecessors, was 
simply an approximation applicable over a limited 
radial range, as noted by King (1978). It is in- 
teresting to note that up until the 1980s, depar- 
tures at large radii from the Reynold's model were 
attributed to tidal-stripping by external gravita- 
tional potentials. That is, for three quarters of a 
century, Reynold's model — originally developed 
from low-quality data for one galaxy — was gener- 
ally thought to describe the original, undisturbed 
stellar distribution in elliptical galaxies. 

de Vaucouleurs' (1948, 1953) i? 1 / 4 surface- 
density model had traction for many years, in 
part due to de Vaucouleurs (1959) arguing that it 
fits better than the Reynold's model used by Hub- 
ble, — a point re-iterated by Kormendy (1977a) 
and others — and the revelation that it fits the 
radially-extended data for NGC 3379 exceed- 
ingly well (de Vaucouleurs & Capaccioli 1979). 
Hodge (1961a, b) had however revealed that dc 
Vaucouleurs' model was inadequate to describe 
faint elliptical galaxies and Hodge (1963, 1964), 
in addition to King (1962 % noted that the 3- 
parameter King model, with its flatter inner pro- 



King (1962) also noted that his model failed to fit the inner 
region of bright elliptical galaxies. 
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file and steeper decline at large radii, did a better 
job. For a time, King's (1962, 1966) model be- 
came popular for describing the light distribution 
in faint elliptical galaxies, at least until the expo- 
nential model — also used for the discs of spiral 
galaxies — was noted to provide a good descrip- 
tion of some dwarf elliptical galaxies (Hodge 1971; 
Faber & Lin 1983; Binggeli et al. 1984) and that 
these galaxies need not have experienced any tidal 
truncation (a prescription of the King model with 
its tidal radius parameter). Lauer (1984, 1985) 
additionally showed that King's modified isother- 
mal model, with its flat inner core, was inadequate 
to describe the deconvolved light-profiles of ordi- 
nary elliptical galaxies with "cores", i.e. galaxies 
whose inner light profile displays a nearly flat core. 
King's model does however remain extremely use- 
ful for studies of star clusters, globular cluster^, 
dwarf spheroidal galaxies and galactic satellites 
which, unlike ordinary elliptical galaxies, can have 
flat cores in their inner surface brightness profile. 

Today, the model of choice for describing 
nearby (and distant) dwarf and ordinary ellipti- 
cal galaxies is Sersic's (1963) generalisation of de 
Vaucouleurs' R 1 / 4 model to give the R x l n surface- 
density model (section I2.1[) . This model repro- 
duces the exponential model when n = 1 and de 
Vaucouleurs' model when n = 4; it can thus de- 
scribe the main body of faint and luminous ellip- 
tical galaxies. The key advantage that this model 
has is (i) its ability to describe the observed stel- 
lar distributions that have a range of central con- 
centrations (known to exist since at least Reaves 
1956) and (ii) it provides a very good description 
of the data over (almost) the entire radial extent. 
Indeed, departures in the light profile from a well- 
fit Sersic's model invariably signal the presence 
of additional features or components, rather than 
any failing of the model. Expanding upon the 
Sersic model, the core-Sersic model (section I2.2j) 
is nowadays used to quantify those galaxies with 
"cores" . 

Although referring to the King model, the fol- 
lowing quote from King (1966) seems particularly 
insightful "... de Vaucouleurs' law appears to refer 
to a particular central concentration and should be 
appropriate only for galaxy profiles that have that 



The Wilson (1975) and Elson (1987) model are also useful 
for describing globular clusters. 



concentration." While noted by others, such as 
Oemler (1976), Capaccioli (1985), Michard (1985) 
and Schombert (1986), some three decades elapsed 
before the relevance of King's remark to ellipti- 
cal galaxies re-surfaced — albeit slowly at first 
- in the 1990s. Indeed, de Vaucouleurs' useful, 
albeit limited, i? 1 / 4 model was referred to as a 
"law" for nearly half a century. However we are 
now more keenly aware that (even normal) ellip- 
tical galaxies possess a range of central concen- 
trations: concentrations which are well quantified 
in terms of the exponent n in Sersic's R}/ n model 
(see Trujillo, Graham & Caon 2001; Graham, Tru- 
jillo & Caon 2001). Although, it should be con- 
fessed that one can still encounter papers which 
use i? 1 / 4 model parameters alongside some model- 
independent measure of galaxy concentration, un- 
aware of the inconsistency arising from the fact 
that every i? 1 / 4 model actually has exactly the 
same level of concentration. 

Before introducing the equation for Sersic's 
model in the following section, it is pointed out 
that in addition to modelling what can be consid- 
ered the main body of the galaxy, one can also 
find excess stellar light at (i) small radii in the 
form of nuclear (i.e. centrally located) discs and 
dense nuclear star clusters (section I2.3[) and also 
at (ii) large radii in the form of halos or envelopes 
in cD and central cluster galaxies (section 12. 4[) . 
As briefly noted above, deficits of stellar flux at 
the cores of massive galaxies are also observed, 
and a model to quantify these stellar distribution, 
relative to the outer non-depleted light profile, is 
described in section 12.21 While non-symmetrical 
components in elliptical galaxies can also exist, 
they are not addressed here given the focus on 
well-structured systems. Somewhat random, non- 
symmetrical components may be a sign of a dis- 
turbed morphology (see E.Barton's Chapter in 
this volume) , of on-going non-uniform star forma- 
tion (see S.Boissier's Chapter in this volume) or 
gravitationally-induced tidal features from exter- 
nal forces. 

2.1. Sersic's model 

Jose Sersic's (1963, 1968) R l / n model, which 
was introduced in Spanish, describes how the pro- 
jected surface-intensity / varies with the projected 
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radius R, such that 
I(R) = I e exp \ -b n 



l/r. 



(1) 



and I e is the intensity at the 'effective' radius R c 
that encloses half of the total light from the model 
(Ciotti 1991; Caon et al. 1993). The term b n 
(« 1.9992n - 0.3271 for 0.5 < n < 10, Capaccioli 
1989) is not a parameter but is instead dependent 
on the third model parameter, n, that describes 
the shape, i.e. the concentration, of the light pro- 
file!!] The exact value of b n is obtained by solving 
the equation T(2n) — 2^{2n, b n ) where 7(2n,x) is 
the incomplete gamma function and Y is the (com- 
plete) gamma function (Ciotti 1991). Useful Sersic 
related expressions have been presented in Ciotti 
(1991), Simonneau & Prada (2004) and Ciotti & 
Bertin (1999), while Graham & Driver (2005) pro- 
vide a detailed review of Sersic's model plus asso- 
ciated quantities and references to pioneers of this 
model. 

The relation between the effective surface 
brightness (fi e = — 2.5 logic) and the central sur- 
face brightness (fio = —2.5 log io) is given by the 
expression 

/i e = Mo + 1-0866, (2) 

where we have dropped the subscript n from the 
term b n for simplicity, while 



(M)c -Mc-2.51og[e f, nr(2n)/6 2 



(3) 



gives the difference between the effective surface 
brightness and the mean effective surface bright- 
ness ((//) c ) within R c . Figure [2] shows the be- 
haviour of the Sersic model. 

Uniting CCD data with wide and deep photo- 
graphic images, Caon et al. (1990, 1993, 1994) 
revealed that the Sersic i? 1 /™ model provided a 
remarkably good description to the stellar distri- 
bution over a large radial range, down to surface 
brightnesses of ~28 B-mag arcsec -2 , for the early- 
type galaxies brighter than Mb = — 18 mag in the 
Virgo cluster. This work was in essence an expan- 
sion of de Vaucouleurs & Capaccioli's (1979) study 
of NGC 3379 which is very well-fit with n = 4. 
Different galaxies were discovered by Caon et al. 



'Ellipticity gradients result in a different Sersic index for the 
major- and minor-axis, as noted by Caon et al. (1993) and 
later quantified by Ferrari et al. (2004). 



(1993) to be equally well fit, but required different 
values of n (see also Bertin et al. 2002). 

Importantly, Caon et al. (1993) additionally 
showed that a correlation existed between stellar 
concentration, i.e. the Sersic index n, and (model- 
independent) galaxy size that was not due to pa- 
rameter coupling in the Sersic model (see also Tru- 
jillo et al. 2001, their section 2). One of the com- 
monly overlooked implications of this result is that 
i? 1 / 4 , and similarly Petrosian (1976), magnitudes, 
sizes and surface brightnesses are systematically in 
error as a function of galaxy concentration (Gra- 
ham et al. 2005; Hill et al. 2011). That is, appli- 
cation of a model which fails to adequately cap- 
ture the range of stellar distributions will result 
in parameters which are systematically biased as 
a function of galaxy mass. For example, fitting an 
i? 1 / 4 model to elliptical galaxies which are actu- 
ally described by an R 1 / 11 model with n less than 
and greater than 4 will yield sizes and luminosi- 
ties which are, respectively, greater than and less 
than the true value (e.g. Binggeli & Cameron 1991; 
Trujillo et al. 2001; Brown et al. 2003). Similarly, 
fitting an exponential model to bulges that are 
best described by an R}l n model with n less than 
and greater than 1 will yield sizes and luminosities 
which are, respectively, greater than and less than 
the true value (e.g. Graham 2001). Obviously one 
does not want to fine tune their galaxy simulations 
to match scaling relations that contain systematic 
biases due to poor measurements, and observers 
are therefore busy fitting R}/ n models these days. 

A good approximation to the internal-density 
profile associated with Sersic's model, i.e. with 
its deprojection, was introduced by Prugniel & 
Simien (1997). Useful expressions for the dynam- 
ics, gravitational potential and forces of this model 
have been developed by Trujillo et al. (2002), 
Tcrzic & Graham (2005) and Terzic & Sprague 
(2007). Somewhat more complex than the early 
light-profile models, such expressions can, impor- 
tantly, accommodate a range of concentrations, 
rather than only varying one scale radius and 
one scale density. Such a model is vital if one 
wishes to properly simulate and understand the 
mass spectrum of elliptical galaxies, whose Sersic 
index n increases with stellar mass. While Gra- 
ham & Driver's (2005) review stated that "No 
attempt has been made here to show the nu- 
merous scientific advances engendered via appli- 
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Fig. 2. — Left panel: Sersic's i? 1 /™ model (equation [T]) for indices n = 0.5,1,2,4 and 10. The effective 
radius (i? e ) and surface brightness (/i c = —2.5 log I e ) have been arbitrarily set to 10 and 25. Right panel: 
Difference between the various surface brightness terms discussed in the text. 



cation of the i? 1 /™ model" , this article reveals how 
Sersic's model, and the core-Sersic model (sub- 
section [22]), have become key in unifying and un- 
derstanding the galaxies around us. 

Like the majority of surface- and internal- 
density models from the last century, the Sersic 
function is an empirical model created to match 
data rather than developed from theory, and as 
such we should be cautious before calling it a 
law. Attempts to find a physical explanation 
for de Vaucouleurs' model yielded results which 
helped to keep it in vogue. Dissipational models 
have long been touted for producing i? 1 / 4 profiles 
(e.g. Larson 1969, 1974), and in the 1980s papers 
based on dissipationless N-body simulations of a 
cold clumpy collapse or the merger of disc galax- 
ies also claimed to finally produce i? 1 / 4 (and also 
Reynold) profiles (e.g. van Albada 1982; McGlynn 
1984; Carlberg et al. 1986; Barnes 1988). How- 
ever a closer inspection reveals clear departures 
from the i? 1 / 4 profile, with the simulated profiles 
better described by an i? 1 /™ model with n < 4. 
Obviously their inability (or perhaps lack of de- 
sire, although see Farouki, Shapiro & Duncan 1983 
whose non-homologous merger remnants were ini- 
tially criticised by i? 1 / 4 aficionados) to create the 
range of stellar concentrations now observed in 
elliptical galaxies highlights a limitation of these 
early works. Nonetheless, these pioneering studies 
have led to N-body simulations by Nipoti et al. 
(2006) and Aceves et al. (2006) — and Farouki et 



al. (1983), whose results with a smaller force soft- 
ening appeared years ahead of their time — have 
now recovered a range of Sersic profile shapes for 
gravitational collapses in a dark matter halo and 
for disc galaxy mergers, respectively. 

Given the empirical nature of Sersic's R 1 / 11 
model, Hjorth & Madsen (1995) revealed how dis- 
sipationless merging and violent relaxation pro- 
vided a physical explanation for the departure 
from the homologous i? 1 / 4 model. Other works 
have explained how the quasi-constant specific en- 
tropy associated with the post violent-relaxation 
stage of elliptical galaxies results in the observed 
mass-dependent range of stellar concentrations in 
elliptical galaxies (Gerbal et al. 1997; Lima Neto 
1999; Marquez et al. 2001). 

It does not seem too unreasonable to spec- 
ulate that elliptical galaxies, whether built by 
near-monolithic collapse, collisions of disc galax- 
ies, wet or dry mergers, appear to eventually ex- 
perience the same force(s) of nature that results 
in their radial stellar distribution depending on 
the total stellar mass. That is, it may not mat- 
ter how the mass was accumulated into an el- 
liptical galaxy, once it becomes a dynamically- 
heated, bound stellar-system, it appears to even- 
tually obey certain universal scaling relations (see 
section [3]) . 

It is interesting to note that Sersic actually 
introduced his model as a way to parameterise 
disc galaxies which he thought were comprised 
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of differing ratios of a disc plus an i? 1 / 4 -bulge. 
His model was not initially intended to fit ellip- 
tical galaxies, and as such it did not immediately 
threaten de Vaucouleurs' model. Credit for popu- 
larising the use of Sersic's i? 1 /™ model for approx- 
imating not only lenticulaiH bulge+disc galaxies 
but for describing pure elliptical galaxies resides 
largely with Massimo Capaccioli (e.g. Capaccioli 
1985, 1987, 1989; Caon et al. 1993; D'Onofrio et 
al. 19940. However, Davies et al. (1988) had also 
introduced this model for dwarf elliptical galax- 
ies, while Sparks (1988) developed an early Gaus- 
sian seeing correction for this model, and Ciotti 
(1991) developed a number of associated expres- 
sions such as the velocity dispersion profile and a 
distribution function. The important quantifica- 
tion that Capaccioli and others provided is how 
the radial distribution of stars in elliptical galax- 
ies, i.e. their concentration, varies with the size, lu- 
minosity and thus the mass of the elliptical galaxy 
(see also Cellone, Forte, & Geisler 1994; Vennik & 
Richter 1994; Young & Currie 1994, 1995; Graham 
et al. 1996; Karachentseva et al. 1996, Vennik et 
al. 1996). As we shall see in this article, the impli- 
cations of this breakthrough have been dramatic, 
unifying what had previously been considered two 
distinct species of galaxy, namely dwarf and ordi- 
nary elliptical galaxies, previously thought to be 
described by an exponential and i? 1 / 4 model, re- 
spectively. 

2.1.1. Dark Matter Halos 

This review would be somewhat incomplete 
without a few words regarding the connection be- 
tween Sersic's model and (simulated) dark matter 
halos. While modified theories of gravity may yet 
make dark matter redundant at some level, it is 
intriguing to note that the Prugniel-Simien (1997) 
internal-density model, developed to approximate 
the deprojected form of Sersic's i? 1 /" model, addi- 
tionally provides a very good representation of the 
internal-density profiles of simulated dark matter 
halos. Merritt et al. (2006) revealed that it actu- 
ally provides a better description than not only the 
Navarro, Frenk & White (1997) model but even a 



9 This term was introduced by Knox Shaw (1915) and 
Reynolds (1920) in their galaxy classification scheme. 

°It is worth noting that D'Onofrio (2001) re-modelled the 
Virgo and Fornax 2-component lenticular galaxies with an 
^/"-bulge plus an exponential disc (see section U 



generalised NFW model with an arbitrary inner 
profile slope 7. 

Sersic's former student, Navarrro, indepen- 
dently applied Sersic's surface-density model to 
the internal-density profiles of simulated dark 
matter halos (Navarro et al. 2004; Merritt et al. 
2005). Jaan Einasto (1965) had previously devel- 
oped this same function as Sersic to describe the 
internal-density profiles of galaxies. Rather than a 
universal profile shape, as advocated by Navarro ct 
al. (1997), a range of simulated dark matter den- 
sity profile shapes is now known to vary with the 
dark matter halo mass (Avila-Reese et al. 1999; 
Jing & Suto 2000; Merritt et al. 2005; Del Popolo 
2010 and references therein). A number of use- 
ful expressions related to this "Einasto model", 
which has the same functional form as Sersic's 
model but is applied to the internal rather than 
projected density profile, can be found in Cardone 
et al. (2005), Mamon & Lokas (2005) and Graham 
et al. (2006). 

An apparent "bulge-halo conspiracy" between 
the radial distribution of stellar mass and dark 
matter (after modification by baryons) has arisen 
in recent years, such that elliptical galaxies re- 
portedly have total internal-density profiles p(r) 
described by power-laws (Bertin & Stiavelli 1993; 
Kochanek 1995; Gavazzi et al. 2007; Buote & 
Humphrey 2011). These power-laws were origi- 
nally claimed to be close to isothermal, such that 
p(r) oc r~ 2 (Koopmans et al. 2006, 2009; Gavazzi 
et al. 2007). Recent work now emphasises that 
only the sample average profile slope is close to 
—2, and that a trend in slope with galaxy size ex- 
ists (Humphrey & Boute 2010; Auger et al. 2010). 
This is a developing field and it is worth not- 
ing that the analyses have been confined to mas- 
sive galaxies with velocity dispersions greater than 
~175 km s _1 , and thus with Sersic indices n > 4 
(Graham, Trujillo & Caon 2001). While the light 
profile shape changes dramatically as the Sersic 
index n increases from ~1 to ~4, there is not such 
an obvious change in light profile shape from ~4 
to higher values of n (see Figure [2]). The apparent 
isothermal profiles of elliptical galaxies, and the 
alleged "bulge-halo conspiracy" , may turn out to 
be a by-product of sample selection (i.e. choosing 
galaxies which have approximately the same struc- 
ture) . It would be interesting to expand the Sloan 
Lens ACS (SLACS) Survey (Bolton et al. 2006) to 
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a greater range than only bright early-type galax- 
ies that are approximately well fit with an i? 1 / 4 
model, and to go beyond the use of simple power- 
laws to describe the total mass density profile once 
the data allows this. Two component mass models 
(Prugniel-Simicn+Einasto) to a range of galaxy 
masses await. Claims that "early- type galaxies are 
structurally close to homologous" may therefore 
be premature, as was the case for the distribution 
of stellar light in elliptical galaxies while the i? 1 / 4 
model was thought to be a law. 

2.2. The core-Sersic model 

The centres of luminous galaxies have long been 
known to possess "cores", such that the surface- 
density profile flattens at the center (e.g. King & 
Minkowski 1966), and King & Minkowski (1972) 
remarked on the inability of the Reynold's and de 
Vaucouleurs' model to match these flattened cores 
in giant galaxies. Although King (1978) identified 
a number of galaxies thought to be well described 
by his model, using seeing-deconvolved ground- 
based images, Lauer (1983, 1984, 1985) analysed 
14 galaxies with 'cores', ranging from 1.5-5.0 arc- 
seconds in radius, revealing that they, like M87 
(Young et al. 1978; Duncan & Wheeler 1980; Bin- 
ney & Mamon 1982), were not exactly described 
by the King model which had a completely flat 
core. Similar conclusions, that cores existed but 
that they do not have flat inner surface brightness 
profiles, were also reported by Kormendy (1982, 
1985a), creating the need for a new model to de- 
scribe the stellar distribution in galaxies. 

Nearly a decade later, the Hubble Space Tele- 
scope was flying and offered factors of a few im- 
provement over the best image resolution achiev- 
able from the ground at that time. Not surpris- 
ingly, astronomers explored the centres of galaxies. 
In an effort to quantify these nuclear regions, af- 
ter the abandonment of the King model and the 
lack of a flattened core in the i? 1 / 4 model, Crane 
et al. (1993), Ferrarese et al. (1994), Forbes et al. 
(1994) and Jaffe et al. (1994) used a double power- 
law model to describe the inner light profiles of 
large galaxies. Grillmair et al. (1994), Kormendy 
et al. (1994) and Lauer et al. (1995) also adopted 
a double-power-law model but one with an addi- 
tional, fifth, parameter to control the sharpness of 
the transition. Their model, which they dubbed 
the "Nukcr law" for describing the nuclear regions 



of galaxies (after excluding any apparent excess 
light), has the same functional form as the double 
power-law model presented by Hernquist (1990, 
his equation 43) to describe the internal-density 
of galaxies (Zhao 1996). 

However, as noted by the above authors, these 
double power-law models were never intended to 
describe the entire radial extent of a galaxy's stel- 
lar distribution, and they provided no connection 
with the outer (i? 1 / 4 -like) radial profile. This dis- 
connection turned out to be their down fall. Due 
to the curved nature of the outer light profiles 
beyond the core, which were being fitted by the 
double power-law model's outer power-law, the 
five parameters of the Nuker model systematically 
changed as the fitted radial extent changed. This 
was first illustrated in a number of diagrams by 
Graham et al. (2003) who revealed that none of 
the Nuker model parameters were robust, and as 
such they could not provide meaningful physical 
quantities. For example, Trujillo et al. (2004) 
reported that the Nuker-derived core-radii were 
typically double, and up to a factor of five times 
larger, than the radius where the inner power-law 
core broke away from the outer i? 1 / 4 -like profile 
- a result reiterated by Dullo & Graham (2011, 
in prep.). An additional problem was that these 
"break radii" were being identified in the so-called 
"power-law" galaxies that showed no evidence of 
a downward departure and flattening from the in- 
ward extrapolation of the outer i? 1 / 4 -like profile. 
This situation arose because of the curved nature 
of what were actually Sersic profiles. That is, 
the so-called "power-law" galaxies not only had 
no distinct "core" like the "core galaxies" do, but 
confusingly they do not even have power-law light 
profiles. 

Given that Caon et al. (1993) and D'Onofrio 
et al. (1994) had established that the Sersic func- 
tion fits the brightness profiles of elliptical galax- 
ies remarkably well over a large dynamic range 
(see the figures in Bertin et al. 2002), it is possi- 
ble to confidently identify departures from these 
profiles that are diagnostic of galaxy formation. 
While in this and the following section we deal 
with partially-depleted cores — also referred to as 
"missing light" — in luminous galaxies (thought 
to be built from dissipationless mergers), the en- 
suing section addresses extra central light above 
the inward extrapolation of the outer Sersic profile 
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(found in galaxies that have experienced dissipa- 
tion and star formation). 

Building on the work of Caon et al. (1993), 
Graham et al. (2003) introduced the core-Sersic 
model, which was applied in Trujillo et al. (2004). 
The model represents a modification to Sersic's 
model such that it has an inner power-law core. 
This represented a dramatic change from what had 
gone before. While the Nuker team (e.g Kormendy 
et al. 1994; Lauer et al. 1995, 2005; Faber et al. 
1997) were combining Nuker model parameters for 
the core with i? 1 / 4 model parameters for the main 
galaxy, Graham et al. (2003), Graham & Guzman 
(2003), Balcells et al. (2003), Graham (2004) and 
Trujillo et al. (2004) advocated the measurement 
of core properties, excesses and deficits of light, 
measured relative to the outer Sersic model. 

For the first time, the core-Sersic model pro- 
vided an expression capable of unifying the nuclear 
regions of galaxies with their outer regions and 
also providing stable physical quantified 1 " 1 !. The 
model can be written as 



I(R) = I' 



1+ ^ 



Rb 

R 



a-i "y/c 



x exp {-b[(R a + i?«)/i? e Q ] 1 /( Q ")}(4) 

where R b denotes the break-radius separating the 
inner power-law having logarithmic slope 7 from 
the outer Sersic function. The intensity I b at the 
break-radius is such that 



I' = I b 2- { ~< /a) exp \b(2 1/a R b /R, 



(5) 



The sixth and final parameter, a, controls the 
sharpness of the transition between the inner 
(power-law) and the outer (Sersic) regime - 
higher values of a indicating a sharper transi- 
tion. Figure [3] shows the core-Sersic model (with 
a = 100) applied to NGC 3348. 

Terzic & Graham (2005) and Terzic & Sprague 
(2007) provide a number of expressions related to 
the potential, force and dynamics^ The core- 



llr The core-Sersic, and also Sersic, model provide robust pa- 
rameters beyond the core if one has sufficient data to sam- 
ple the curvature in the light profile, in practice this re- 
quires radial data out to ~li? e . 

12 The complex nature of this model has resulted in the ap- 
pearance of alternate expressions, e.g. Spergel (2010, see 
his figure 3). 



Sersic model is further discussed and used by Fer- 
rarese et al. (2006a,b), Cote et al. (2006, 2007), 
Kawata, Cen & Ho (2007) and Ciotti (2009). 

2.2.1. Central Mass Deficits 

The collisional construction of galaxies from the 
merger of lesser galaxies is thought to be a com- 
mon occurrence in the Universe. Coupled with 
the presence of a supermassive black hole (SMBH) 
at the heart of most galaxies (Wolfe & Burbidge 
1970; Magorrian et al. 1998), dissipationless merg- 
ers were proposed by Begelman, Blandford & Rees 
(1980; see also Ebisuzaki, Makino, & Okumura 
1991) to explain the depleted nuclei, i.e. the cores, 
observed in giant elliptical galaxies (e.g. King 
1978, and references therein). It is thought that 
core-depletion is primarily due to the gravitational 
slingshot (Saslaw, Valtonen, & Aarseth 1974) ef- 
fect that the coalescing SMBHs — from the pre- 
merged galaxies — have on stars while they them- 
selves sink to the bottom of the potential well of 
the newly wed galaxvP^I 

Theory predicts that the central mass deficit 
Md e f should scale with 0.5iVMbh, where Afbh is 
the final (merged) black hole mass and N the num- 
ber of major "dry" (i.e. gas free, dissipationless) 
mergers (Milosavljevic & Merritt 2001; Merritt 
k Milosavljevic 2005; Merritt 2006a,b). Graham 
(2004) used the core-Sersic model to quantify the 
central deficit of stars relative to the inward ex- 
trapolation of the outer Sersic profile. Figure [4] 
suggests that the luminous elliptical galaxies sam- 
pled have experienced an average of 1 or 2 major 
dry (i.e. dissipationless) mergers, a conclusion in 
agreement with select ACDM models of galaxy for- 
mation (Haehnelt & Kauffmann 2002; Volonteri et 
al. 2003). 

Quantification of the central stellar deficit rela- 
tive to the inward extrapolation of the outer Sersic 
profile has also been applied to bright Virgo clus- 



13 The presence of FUV—NUV colour gradients in core galax- 
ies such as NGC 1399 suggests that they may not have 
been built from major, dry merger events (Carter et al. 
2011). Additionally, the globular cluster specific frequency 
in core galaxies may be at odds with core-galaxy forma- 
tion through equal-mass merger events because the fainter, 
intermediate luminosity elliptical galaxies may have lower 
specific frequencies (Harris & van den Bergh 1981). Alter- 
native ideas for core-formation in giant galaxies have been 
proposed: Boylan-Kolchin et al. 2004; Nipoti et al. 2006; 
Martizzi et al. 2011. 
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Fig. 3. — Left: Core-Sersic model (solid line) fit to the major-axis i?-band light profile of NGC 3348 (dots), 
with the dashed line showing the associated Sersic component. The inner depleted zone corresponds to a 
stellar mass deficit of 3 x 10 8 Mq (see Graham 2004). Right: NGC 5831 has, in contrast, no (obvious) 
partially depleted core and is well described by the Sersic's model alone. The rms scatter is shown in the 
lower panels. Figure taken from Graham et al. (2003). 
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Fig. 4. — Central mass deficit for seven core galax- 
ies derived using the core-Sersic model (circles) 
and the Nuker model (stars) plotted against each 
galaxy's predicted central supermassive black hole 
mass. The solid and dashed line shows Md e f equal 
1 and 2 Mbh, respectively. Figure adapted from 
Graham (2004). 



ter galaxies by Ferrarese et al. (2006a) and Cote et 
al. (2007) and with the exception of VCC 798 — a 
lenticular (bulge plus disc) galaxy — provides sim- 
ilar results. Of course when an outer disc exists, a 
core-Sersic bulge plus disc fit is required, otherwise 
the disc will bias the Sersic parameters P^l 

2.3. Excess nuclear light 

While galaxies brighter than Mb ~ —20.5 mag 
have partially depleted cores (e.g. Faber et al. 
1997; Rest et al. 2001; Graham & Guzman 2003), 
fainter galaxies often contain additional stellar 
components at their centres. Such excess nuclear 
light, above that of the underlying host galaxy, 
has long been known to exist (e.g. Smith 1935; 
Reaves 1956, 1977; Romanishin, Strom & Strom 
1977) and was systematically studied by Binggeli 



Trujillo ct al. (2004) and Graham (2004) intentionally ex- 
cluded lenticular galaxies from their publication because 
it was felt that the community were not yet ready for an 
8-parameter model (6 core-Sersic plus 2 exponential disc 
parameters). Intriguingly, VCC 798 resides near the 1-to-l 
line in Figure[4]when an exponential-disc plus (core-Sersic)- 
bulge fit is performed (Graham, unpublished work). 
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et al. (1984), van den Bergh (1986) and Binggeli 
& Cameron (1991, 1993) in a number of dwarf 
elliptical galaxies. As far back as Larson (1975) 
it was known that simulations containing gas can 
account for these dense nuclear star clusters; clus- 
ters which became easier to detect with the Hub- 
ble Space Telescope (e.g. Carollo, Stiavelli & Mack 
1998). 

For many years it was common practice to sim- 
ply exclude these additional nuclear components 
from the analysis of the galaxy light profile (e.g. 
Lauer et al. 1995; Byun et al. 1996; Rest et al. 
2001; Ravindranath et al. 2001). Using HST data, 
Graham & Guzman (2003) simultaneously mod- 
elled the host galaxy and the additional nuclear 
component with the combination of a Sersic func- 
tion for the host galaxy plus a Gaussian func- 
tion for the nuclear star cluster. As we will see 
in section [H they also showed that the lenticular 
galaxies in their early-type galaxy sample could be 
modelled via a Sersic-bulge plus an exponential- 
disc plus a Gaussian- (star cluster) decomposition 
of their light profile — as done by Wadadekar et 
al. (1999) with ground-based images. Many other 
studies have since modeled the nuclear star clus- 
ters seen in HST images, see Figure [5J with the 
combination of a nuclear component plus a Sersic 
host galaxy (e.g. Grant et al. 2005; Cote et al. 
2006; Ferrarese et al. 2006a; Graham & Spitler 
2009). 

While Graham & Guzman (2003) and Cote et 
al. (2006) found that some nuclear star clusters 
could actually be resolved, it is not yet established 
what mathematical function best describes their 
radial distribution of stars. The closest example 
we have to study is of course the 30 million so- 
lar mass nuclear star cluster at the centre of the 
Milky Way (Launhardt et al. 2002). Graham & 
Spitler (2009) provided what may have been the 
first ever analysis after allowing for the significant 
contamination from bulge stars. Although they 
found that the cluster could be well described by 
a Sersic index with n = 3 (see Figure 12 .31 it re- 
mains to be tested how well a King model can 
describe the data, or at least the old stellar pop- 
ulation known to have a core (e.g. Genzel et al. 
1996). The excess nuclear light in M32 has also 
been well fit with an n — 2.3 Sersic function (Gra- 
ham & Spitler) but this too may yet be better 
described by a King model. What is apparent is 



that neither the underlying bulge nor the nuclear 
component are described by a power-law. Theo- 
ries which form or assume such power-law cusps 
appear to be at odds with current observations. 

2.4. Excess halo light 

This section shall be somewhat cursory given 
R. Bower's detailed chapter on clusters and the in- 
tracluster medium in this volume. 

Brightest cluster galaxies (BCGs), residing 
close to or at the centres of large galaxy clus- 
ters, have long been recognised as different from 
less luminous elliptical galaxies: their light pro- 
files appear to have excess flux at large radii (e.g. 
Oemler 1976; Carter 1977; van den Bergh 1977; 
Lugger 1984; Schombert et al. 1986). However, 
before exploring this phenomenon, it is important 
to recall that the light profiles of large galaxies 
have Sersic indices n greater than 4 (e.g. Caon et 
al. 1993; Graham et al. 1996)0 Subsequently, 
at large radii in big elliptical galaxies, there will 
be excess flux above that of an i? 1 / 4 model fit to 
some limited inner radial range. 

An initially puzzling result from the Sloan Digi- 
tal Sky Survey (SDSS; York et al. 2000) survey was 
the lack of light profiles with Sersic n > 5-6 (Blan- 
ton et al. 2005a). This was however soon resolved 
when Blanton et al. (2005b), Mandelbaum et al. 
(2005) and Lisker (2005, 2006b, 2007) pointed out 
a serious sky-subtraction problem with the early 
SDSS Photometric Pipeline (photo: Ivezic et al. 
2004). The sky- value to be subtracted from each 
galaxy had been measured to close to the galaxy 
in question, and because galaxies with large Sersic 
indices possess rather extended light-profiles, their 
outer galaxy light was actually measured and then 
subtracted from these galaxies. As a result, the 
high Sersic index light-profiles were erroneously 
erased and missed from the SDSS. This resulted 
in the i? 1 / 4 model appearing to provide good fits 
to bright elliptical galaxies, and consequently a se- 
ries of papers based on structurally biased radii, 
magnitudes and surface brightnesses. 

Bearing in mind that large elliptical galaxies 
have high Sersic indices, it is important to dis- 
tinguish between (i) an inadequacy of the i? 1 / 4 
model to describe what is actually a single R}/ n 



5 As n — > oo, the Sersic model can be approximated by a 
power-law (see Graham & Driver 2005). 
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Fig. 5. — Left: 5-parameter Nuker model fit to the y-band light-profile of the nucleated galaxy NGC 596 after 
excluding the inner three data points (Lauer et al. 2005). Right: 3-parameter Sersic model plus 2-parameter 
point-source (Gaussian) fit to the same light profile of NGC 596. With the same number of parameters, 
this model fits both the inner, intermediate, and outer light-profile. Figure from Dullo & Graham (2011, in 
prep.). 
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Fig. 6.— Uncalibrated, 2MASS, J^-band inten- 
sity profile from the centre of the Milky Way, taken 
from Schodel et al. (2009, their Figure 2). The 
nuclear star cluster is modelled with a Sersic func- 
tion with n — 3 (dotted curve) and the underlying 
host bulge with an exponential function (Kent et 
al. 1991) that has an effective half-light radius of 
~4.5 degrees (e.g. Graham & Driver 2007) and 
is therefore basically a horizontal line. One par- 
sec equals 25 arcseconds. Figure from Graham & 
Spitler (2009). 



profile, and (ii) a distinct physical component such 
as an envelope of diffuse halo light surrounding a 
central galaxy. Early quantitative photometry of 
cD galaxies (supergiant D galaxies: e.g. Matthews 
et al. 1964; Morgan & Lesh 1965) suggested the 
presence of an inner i? 1 / 4 spheroid plus an outer 
exponential corona (de Vaucouleurs 1969; de Vau- 
couleurs & de Vaucouleurs 1970). One should 
however question if this outer corona is a distinct 
entity or not. To answer this in the affirmative, 
astronomers can point to how the light profiles 
can display inflections marking the transition from 
BCG light to intracluster light. Gonzalez et al. 
(2005) additionally showed that the inflection in 
the light profiles of many BCGs was also associ- 
ated with an a change in the ellipticity profile, sig- 
nalling the switch from BCG light to intracluster 
light. 

Gonzalez et al. (2005) and Zibetti et al. (2005) 
chose to model their BCG sample using an i? 1 / 4 + 
i? 1 / 4 model to describe the inner galaxy plus the 
outer halo light. However, given that elliptical 
galaxies are better described by the R}l n model, 
and the desire to measure the actual concentra- 
tion of halos rather than assign a fixed i? 1 / 4 profile, 
Seigar et al. (2007) fitted an R 1 ' n + R 1 / n model to 
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the light profiles of five BCGs. Not surprisingly, 
they found that the R 1/4 + R 1/4 model was not 
the optimal representation. An i? 1 /™-galaxy plus 
an exponential-halo model was found to provide 
the optimal fit in three instances, with an addi- 
tional galaxy having no halo detected. The associ- 
ated galaxy-to-halo luminosity ratios can be found 
there. This rc-rcvclation of an exponential model, 
rather than an i? 1 / 4 model, describing the halo has 
since been confirmed by Pierini et al. (2008). In- 
triguing is that the halo does not trace the NFW- 
like dark-matter halo density profiles produced in 
ACDM simulations (section 12.1. 1[) . Stellar halos 
around non-BCG galaxies have also now been re- 
ported to display an exponential radial distribu- 
tion (e.g. Gadotti 2011; Tal & van Dokkum 2011). 

3. Structure related scaling relations 

While it is common practice, and somewhat 
helpful, to call elliptical galaxies fainter than 
about Mb = —18 mag by the term "dwarf el- 
liptical" (Sandage & Binggeli 1984), it should be 
noted, and this section will reveal, that on all 
measures they appear to be the low-mass end of 
a continuous sequence which unifies dwarf and 
normal/luminous elliptical galaxies. Not only is 
this true in terms of their structural properties 
(e.g. Binggeli et al. 1984; Graham & Guzman 
2003; Gavazzi et al. 2005; Ferrarese et al. 2006a; 
Cote et al. 2006, 2007, 2008; Misgeld et al. 2008, 
2009; Janz & Lisker 2009; Graham 2010; Chen et 
al. 2010; Glass et al. 2011) but even the degree 
of kinematically distinct components is similar 
(Chilingarian 2009). 

There are many important relations between 
stellar luminosity, colour, metallicity, age, and dy- 
namics that reveal a continuous and linear be- 
haviour uniting dwarf and giant elliptical galaxies 
(e.g. Caldwell 1983; Davies et al. 1983; Binggeli 
et al. 1984; Bothun et al. 1986; Geha et al. 2003; 
Lisker & Han 2008). However, Kormendy et al. 
(2009, their section 8) dismiss all of these appar- 
ently unifying relations by claiming that they must 
not be sensitive to different physical processes 
which they believe have produced a dichotomy be- 
tween dwarf and ordinary elliptical galaxies. They 
claim that it is only relations which show an appar- 
ent different behavior at the faint and bright end 
that are sensitive to the formation physics and the 



remainder are not relevant. This section explains 
why such non-linear relations are actually a conse- 
quence of the (dismissed) linear relations, and as 
such these non-linear relations actually support a 
continuum between dwarf and ordinary elliptical 
galaxies. 

To begin, it should be reiterated that (dwarf 
and ordinary) elliptical galaxies — and the bulges 
of disc galaxies fsection [4.1l) — do not have struc- 
tural homology. Instead, they have a continuous 
range of stellar concentrations — quantified by 
the Sersic index n (Davies et al. 1988; Caon et 
al. 1993; D'Onofrio et al. 1994; Young & Curric 
1994, 1995; Andredakis et al. 1995) — that varies 
linearly with both stellar mass and central surface 
brightness (after correcting for central deficits or 
excess light). A frequently unappreciated conse- 
quence of these two linear relations is that rela- 
tions involving either the effective half-light radius 
(i? e ) or the effective surface brightness (/i c ), or the 
moan surface brightness within R c ((/i) c ), will be 
non-linear. Such curved relations have often been 
heralded as evidence that two different physical 
processes must be operating because the relation 
is not linear and has a different slope at either end. 
To further complicate matters, sample selection 
which includes faint and bright elliptical galax- 
ies, but excludes the intermediate-luminosity pop- 
ulation, can effectively break such continuously 
curved relations into two apparently disconnected 
relations, as can selective colour-coding. 

There are three distinct types of (two-parameter) 
relations involving the properties of elliptical 
galaxies: (i) linear relations which are taken to re- 
veal the unified nature of dEs and Es; (ii) curved 
relations revealing a continuity that had in the 
past been mis-interpreted to indicate that distinct 
formation process must be operating; and (iii) 
broken relations which imply that two physical 
mechanisms are operating. In the following sec- 
tions we shall learn how the linear relations result 
in the existence of curved relations when using 
effective radii and surface brightnesses. We shall 
also see that the transition in the broken relations 
occurs at Mb ~ —20.5 mag and thus has nothing 
to do with the previously held belief that dEs and 
Es are two distinct species separated at Mb = — 18 
mag (Wirth & Gallager 1984; Kormendy 1985b, 
2009). 
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3.1. Linear relations 

This section introduces two key relations, from 
which a third can be derived, involving struc- 
tural parameters. They are the luminosity- 
concentration (L-ri) relation and the luminosity- 
(central density j3 (L-u ) relation. We shall use 
the Sersic shape parameter n to quantify the cen- 
tral concentration of the radial light profile, and 
use the projected central surface brightness /iq as 
a measure of the central density^ 

It is noted that one can expect a certain level 
of scatter in the L-n and L-fio diagrams because 
both the central density and the radial concentra- 
tion of stars that one observes depends upon the 
chance orientation of one's triaxial galaxy (Bin- 
ney 1978). This is of course also true for mea- 
surements of effective surface brightness, half-light 
radii, velocity dispersions, etc. To have the best 
relations, it is important that we use Sersic param- 
eters from elliptical galaxies rather than parame- 
ters from single-Sersic fits to samples of elliptical 
and lenticular galaxies given the two-component 
(2D-disc plus 3D-bulge) nature of the lenticular 
galaxies. Given the offset nature of bulges and 
elliptical galaxies in the L-n diagram (e.g. Gra- 
ham 2001; see also Mollenhoff & Heidt 2001) it 
is also important that bulges not be combined in 
this section's analysis of elliptical galaxy scaling 
relations. 

3.1.1. Luminosity- (central surface brightness) 
relation 

Caldwell (1983; his Figure 6) and Bothun et al. 
(1986, their figure 7) revealed that, fainter than 
Mb ~ —20.5 mag, there is a continuous linear 
relation between luminosity and central surface 
brightness. Furthermore, Binggeli et al. (1984, 
their figure 11) and Binggeli & Cameron (1991, 
their figure 9 and 18) revealed that, when using 
the inward extrapolation of King models, the L— 
Ho relation is continuous and roughly linear from 
— 12 > Mb > —23 mag. This same general result 
was also highlighted by Jerjen & Binggeli (1997) 
and Graham & Guzman (2003) when using the in- 



°Here the "central density" refers to the density prior to 
core depletion in giant elliptical galaxies or the growth of 
additional nuclear components in smaller elliptical galaxies. 
17 To convert from the surface density to the internal density, 
one can use equation 4 from Terzic & Graham (2005). 



ward extrapolation of the Sersic model. The ben- 
efit of this approach is that one's central surface 
brightness is not biased by the presence of a de- 
pleted core or any additional nuclear components 
within the host galaxy. Figure [T^i displays the el- 
liptical galaxy (Mb, Ho) data set from Graham & 
Guzman (2003) fit by the expression 

M B = 0.67m , b - 29.5. (6) 

The actual central surface brightness of the 
luminous "core galaxies" is shown in Figure [7^, 
rather than the value obtained from the inward 
extrapolation of their outer Sersic profile. As such 
these "core galaxies" were excluded from the fit, 
but see the discussion in section [3~3l As an aside, 
if the central supermassive black hole mass -Mbh 
in elliptical galaxies is directly related to the cen- 
tral stellar density (see Graham & Driver 2007), 
then the connections between Mbh and the global 
galaxy properties, such as total mass and velocity 
dispersion, may be secondary. 

3.1.2. Luminosity- concentration relation 

The linear relation between luminosity and 
Sersic index, or strictly speaking the logarithm 
of these quantities, has been shown many times 
(e.g. Young & Currie 1994; Jerjen & Binggeli 1997; 
Graham & Guzman 2003; Ferrarese et al. 2006a). 
This continuous relation between magnitude and 
concentration^! for elliptical galaxies had of course 
been recognised before (e.g. Ichikawa et al. 1986, 
their figure 11). The following Ms~n expression 
is shown in Figure |7fc, again matched to the sam- 
ple of elliptical galaxies compiled by Graham & 
Guzman (2003). 

M B = -9.4 log(n) - 14.3. (7) 

Graham & Guzman (2003) excluded two-component 
lenticular galaxies fit by others with a single- 
component Sersic model. It may be prudent to 
continue to exclude these galaxies even after a 
Sersic-bulge plus exponential disc fit because the 
Ms~n relation defined by bulges, at least in spiral 
galaxies, is different to that defined by elliptical 
galaxies (Graham 2001, his figure 14). 



Graham et al. (2001) contains a review of various concen- 
tration indices used over the decades, while Trujillo et al. 
(2001) was the first to quantify the monotonic relation be- 
tween Sersic index and concentration. 
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Fig. 7. — Linear relations between the observed B-band central surface brightness (ho,b) with a) the absolute 
B-band magnitude Mb and b) the Sersic index n for a sample of elliptical galaxies. Panel c) shows the galaxy 
magnitudes versus the Sersic indices, with the line given by equation [7] The "core galaxies" (large filled 
circles) with partially depleted cores can be seen to have lower central surface brightnesses than the relation 
in panel a) . Inward extrapolation of their outer Sersic profile yields fiQ values which follow the linear relation, 
as previously noted by Jerjen & Binggeli (1997). The data have come from the compilation by Graham & 
Guzman (2003, their figure 9). Dots represent dE galaxies from Binggeli & Jerjen (1998), triangles represent 
dE galaxies from Stiavelli et al. (2001), large stars represent Graham & Guzman's (2003) Coma dE galaxies, 
asterix represent intermediate to bright E galaxies from Caon et al. (1993) and D'Onofrio et al. (1994), open 
circles represent the so-called "power-law" E galaxies from Faber et al. (1997) and the filled circles represent 
the "core" E galaxies from these same Authors. 



3.1.3. Concentration- (central surface brightness) 
relation 

Combining the above two equations provides an 
expression between central surface brightness and 
Sersic index such that 



Ha = 22.8 - 14.11og(n), 



(8) 



which is shown in Figure [7jj, where it can be seen 
to be roughly applicable for values of n > 1 . 

3.2. Curved relations 

This section explains why, as a direct result of 
the above linear relations — which unite dwarf 
and giant elliptical galaxies — expressions involv- 
ing either the effective half-light radius i? e , the as- 
sociated effective surface brightness [i c at this ra- 
dius, or the mean surface brightness (/Lt) e enclosed 
within this radius, are curved. 

3.2.1. Luminosity- (effective surface brightness) 
relation 

The following analysis is from Graham & 
Guzman (2003). 



Given the empirical Mg-n relation (equa- 
tion [7]), one knows what the expected value of 
n is for some value of Mb- One can then convert 
the empirical Ms-fJ-o relation (equation [5]) into an 
Ms-f^e relation using the exact relation between 
IIq and /Li e which depends only on the value of n 
(equation^. Doing so, one obtains the expression 

fi c = 1.5M B + 44.25 +1.0866, 

= 1.5[M B + 14.3] +22.8 + 1.0866, (9) 

where b as 1.9992rt — 0.3271 and equation [7] is 
used to replace n in terms of Mb, such that 

n = 10 -(14.3+M B )/9.4_ 

One can similarly convert the empirical Mb~ 
relation (equation [S]) into an Ms-(/i) e relation 
using the exact relation between /io and (/i) c which 
also depends only on the value of n (equation [3]). 
Doing this, one obtains the expression 



(/i) c = 1.5A/ B +44.25+1.0866-2.51og 



e b nT(2n) 



b 2n 



(10) 

where again b sa 1.9992n — 0.3271 and equation [7] 
is used to replace n in terms of Mb ■ These curves 



17 



can be seen in Figure [8] (adapted from Graham & 
Guzman 2004). 

Binggeli et al. (1984, their figure 8) and Capac- 
cioli & Caon (1991) previously showed with em- 
pirical data that the Mg-(/i) e relation is curved. 
What was new in Graham & Guzman (2003) was 
the explanation. In the past, evidence of non- 
linear relations involving parameters from dwarf 
and giant elliptical galaxies were heralded as evi- 
dence of a dichotomy. Coupled with galaxy sam- 
ple selection that excluded the intermediate pop- 
ulation, and therefore resulted in two apparently 
disconnected relations, acted to further convince 
some that they were dealing with two classes of 
ObjectPal 

3.2.2. Size- Luminosity relation 

Now that we know how to play this game, 
one can additionally make predictions for rela- 
tions involving the effective radius R e because we 
know that the luminosity L = 2n(I) e R%,, with 
(ju) e = — 2.51og(J) e - As explained in Graham & 
Worley (2008, their section 5.3.1), one can derive 
the size-luminosity relation such that 

logi? e [kpc] = — f + 1.066+0.434n+0.51og - 1 — - 
10 [e°nl (^ 

(11) 

for 0.5 < n < 10, with b « 1.9992n - 0.3271 and 
where equation [7] is used to replace n in terms 
of Mb- This size-luminosity relation for ellipti- 
cal galaxies is shown in Figure [5J: along with real 
galaxy data. 

Binggeli et al. (1984, their figure 7; cf. Misgeld 
& Hilker 2011, their figure 1) also demonstrated, 
with empirical data, that the L — R e relation for 
dwarf and giant elliptical galaxies is curved. Their 
diagram, in addition to Figure [5J; seen here, re- 
veals why studies which only sample bright ellip- 
tical galaxies are often contempt to simply fit a 
straight line (e.g. Kormendy 1977b). The expla- 
nation for why this happens is of course akin to 
approximating the Earth as flat when one (forgiv- 
ably) does not sample enough of what is actually 



! A simple one-dimensional example of sample bias would be 
a survey of the average physical properties, such as size or 
mass, of people at a primary school. One would measure 
the properties of children and adults, but miss the bridging 
population which reveals a continuity and thus unification 
of the species. 



a curved profile. As Graham et al. (2006, their fig- 
ure lb) re-revealed recently, and as reiterated by 
Bernardi et al. (2007), a sample of massive ellipti- 
cal galaxies will have a steeper size-luminosity re- 
lation than a sample of ordinary elliptical galaxies, 
which will in turn have a steeper size-luminosity 
relation than a sample of dwarf elliptical galax- 
ies because the size-luminosity relation is curved. 
Graham & Worley (2008) explains why the L — R e 
relation given by equation (TTJ based on two linear 
relations and the functional form of Sersic's model, 
is curved. 

Interestingly, due to the linear relations be- 
tween magnitude, central surface brightness and 
the logarithm of the Sersic exponent n, the use of 
faint isophotal radii results in what is roughly a 
linear size-luminosity relation (e.g. Oemler 1976; 
Strom & Strom 1978; Forbes et al. 2008; van den 
Bergh 2008; Nair et al. 2011), with the bright-end 
slope dependent on the adopted isophotal limit or 
Petrosian radius used. The implications of this 
important observation shall be detailed elsewhere. 

Helping to propagate the belief that dwarf and 
ordinary elliptical galaxies are distinct species, 
Dabringhausen et al. (2008) and Lisker (2009) fit 
a double power-law to their curved size-luminosity 
relation for dwarf and ordinary elliptical galaxies, 
thus yielding distinct slopes at the faint and bright 
end. In addition, the interesting study by Janz & 
Lisker (2008) reported small deviations from the 
predicted curved relation. However, galaxies that 
are well described by Sersic's function and which 
follow linear M—n and M-/xo relations must fol- 
low a single curved M-R e relation. The devia- 
tions that they found are therefore mirroring a) 
the inadequacy of the fitted linear relations to the 
M—n and M-/iq distribution (a point noted in the 
caption of Figure |8]) and/or b) poor fitting Sersic 
models to their sample of elliptical and disc galax- 
ies. Adding uncertainties to the linear relations in 
section I3H1 and propagating those through to the 
predicted M-R c relation is required before we can 
claim evidence of significant deviations. However, 
searching for such second order effects may indeed 
be interesting given the different types of dwarf 
galaxies that are emerging (Lisker et al. 2007). 

3.2.3. Size- concentration relation 

One can additionally derive an expression relat- 
ing R c and n by substituting the magnitude from 



18 



-22 
r w -20 

US 

J -16 
-14 
-12 



1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 
a) L^... J 


_i i i i i i i i i i i i i i i i i i i_ 
Lb) i 


_ 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 _ 


- 


- 


— ° **>^^*^ °x — 

- i .- 



28 26 24 22 20 

<M>e 



26 24 22 20 
Me 



-0.5 0.5 1 
R [kpc] 



1.5 



Fig. 8. — Elliptical galaxy B-band magnitude versus a) mean effective surface brightness (equation ITO)) . 
b) effective surface brightness (equation [9J , and c) effective radius (equation I lip. The continuous, curved 
relations are predictions from the observed linear relations in Figure [71 The slight mismatch seen here arises 
from the rough fitting-by-eye of the linear relations. The data points have come from the compilation by 
Graham & Guzman (2003, their figure 9). 



the empirical M^-n relation, expressed in terms 
of n (equation [7]), into the size-luminosity relation 
(equation [TTj) to give 



log R e [kpc] = 



0.434n- 0.364- 0.94 log(n 

b 2n 



-0.5 loe 



e b nT(2n) 



(12) 



While the Ms~n relation is linear, the R e -n rela- 
tion is curved, as can be seen in Figure [S] 

In passing it is noted that the form of this re- 
lation (equation [12]) matches the bulge data from 
Fisher k Drory (2010, their Figure 13). They in- 
terpret the departure of the low-n bulges (n < 2) 
from the approximately linear relation defined by 
the high-n bulges (n > 2) to indicate that a differ- 
ent formation process is operating to produce the 
less concentrated "pseudobulges" . However, based 
upon linear unifying relations that span the artifi- 
cial n = 2 divide, we know that this R e -n relation 
must be curved. Without an understanding of this 
relation, and other curved relations (e.g. Greene, 
Ho & Barth 2008), they have at times been misin- 
terpretted and used to claim the existence of dif- 
ferent physical processes (see section I4T31 for a dis- 
cussion of pseudobulges: Hohl 1975 and references 
therein). 

It may be worth better defining the behavior of 
the R e -n relation at small sizes in Figure [9l The 
data from Davies et al. (1988) suggests than when 
n = 0.5, values of R c may range from 1 kpc down 



to 0.2 kpc (Caon et al. 1993, their figure 5). Such 
a reduction to the flattening of the R c -n distri- 
bution, below n « 1, may in part arise from the 
inclusion of dwarf spheroidal galaxies (see Misgeld 
& Hilker 2011, their figure 1). 

3.2.4- Size- (effective surface brightness) relation 

As discussed in Graham (2010), the first two 
linear relations in Figure [7] naturally explain the 
curved (ii) e —R e relation in Figure [T0l From the 
empirical /io -n relation (equation (S[ Figure [TJd) , 
one can convert /io into (//) e using equations [5] and 
[3] The effective radius R c is acquired by match- 
ing the empirical Mb— Ho relation (equation[6j Fig- 
ure [7^) with the absolute magnitude formula 

M = (n) B - 2.5 log(27ri? 2 kpc ) - 36.57, (13) 

(see Graham & Driver 2005, their equation 12). 
Eliminating the absolute magnitude gives the ex- 
pression 



logi? c 



1 f (l«>e 



9.07 + 0.726- 1.67 log 



nc b T(2n) 
(14) 



in which we already know the value of n associ- 
ated with each value of {fx) e - This is achieved by 
(again) using the empirical fJ,o~n- relation (equa- 
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tion [8]) with equations [2] and |3l such that 



(fi) c = 22.8+1.0866-14.1 log(n)-2.5 log 



neT(2n) 



(15) 

and 6 » 1.9992n - 0.3271. Equation O obtained 
from two linear relations involving Sersic parame- 
ters, is a curved relation that is shown in Figure lTOl 
Overplotted this predicted relation are data points 
from real galaxies. 

For those who may have the Sersic parameter 
set (i£ e ,// e ,n), one can use equation [3] to convert 
fi c into (/i) if one wishes to compare with the 
relation given by equation 1 141 For those who may 
have the parameter set (R e , // e ), perhaps obtained 
with no recourse to the Sersic model, equation 1141 
can easily be adjusted using equation [3] to give a 
relation between R c and fi c such that 



logi? c 



Me 

15 



1.81 -0.5 log 



ne b T(2n) 

pgO.296 fj2n 



(16) 



where the value of n associated with the value of 
fi c is given by 

/i e = 22.8 + 1.0866- 14.1 log(n). (17) 

To summarise, due to the linear relations in Fig- 
ure [7] which connect dwarf and ordinary elliptical 
galaxies across the alleged divide at Mb = —18 
mag (Kormendy 1985), or at n = 2 (Kormendy & 
Kennicutt 2004), coupled with the smoothly vary- 
ing change in light profile shape as a function of 
absolute magnitude, the (/z) e --R e and /J e -Re rela- 
tions are expected to be curved (FigureHU]). as pre- 
viously shown with empirical data by, for example, 
Capaccioli & Caon (1991). This also explains why 
the fitting of a linear relation to (R c , (J-c) data by 
Hoessel et al. (1987) resulted in slopes that de- 
pended on their galaxy sample magnitude. 

The Kormendy relation is a tangent to the 
bright arm of what is actually a curved distribu- 
tion defined by the relation given by equation [T4l 
that is taken from Graham (2010). The appar- 
ent deviant nature of the dwarf elliptical galax- 
ies from the approximately linear section of the 
bright-end of the (fj) e -R c distribution does not ne- 
cessitate that a different physical process be oper- 
ating. Moreover, as noted by Graham & Guzman 
(2004) and Graham (2005), galaxies which appear 
to branch off at the faint end of the Fundamen- 
tal Plane (Djorgovski & Davis 1987) — the flat 
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The dashed line shows the Mb = —13 mag limit 
for the Virgo cluster. 
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portion at the bright end of a curved hypersurface 
— also need not have formed from different phys- 
ical mechanisms. Simulations that assume or re- 
produce a linear (fj) e -R e or n e -R relation, across 
too great a magnitude range, have failed to mimic 
the continuous curved distribution defined by real 
elliptical galaxies. The same remark is true for 
simulations of the 'Fundamental Plane'. 

3.3. Broken relations 

3.3.1. Luminosity- (central surface brightness) 
relation 

While the relation between a galaxy's absolute 
magnitude and extrapolated central surface bright- 
ness is remarkably linear (section f3.1.1|) . there is a 
clear break in this relation when using the actual 
central surface brightness at the luminous end of 
the distribution (Figure [7h). This departure from 
the Mb-^o relation by elliptical galaxies brighter 
than M B « -20.5 mag (M > 0.5-1 x 10 n M Q ) 
was addressed by Graham & Guzman (2003) in 
terms of partially depleted cores relative to the 
outer Sersic profile (see also Graham 2004; Tru- 
jillo et al. 2004; Merritt k Milosavljevic 2005; 
Ferrarese et al. 2006a; Cote et al. 2007). This 
transition has nothing to do with the alleged di- 
vide between dwarf and giant elliptical galaxies at 
around Mb = — 18 mag, but is instead directly re- 
lated with the Sersic versus core-Sersic transition 
at around Mb — —20.5 mag. 

As noted in section I2.2.1) such partially de- 
pleted cores in luminous core-Sersic galaxies are 
thought to have formed from dry, dissipationless 
galaxy merger events involving the central coa- 
lescence of supermassive black holes (Begelman, 
Blandford, & Rees 1980; Ebisuzaki, Makino, & 
Okumura 1991; but see footnote fT3|) and resulted 
in Trujillo et al. (2004) advocating a "new ellipti- 
cal galaxy paradigm" based on the presence of a 
central stellar deficit versus either none or an ex- 
cess of light, an approach embraced by Ferrarese 
et al. (2006a), Cote et al. (2007) and others. 

Further evidence for a division at Mb = —20.5 
mag comes from the tendency for the brighter 
galaxies to be anisotropic, pressure supported el- 
liptical galaxies having boxy isophotes, while the 
less luminous early-type galaxies may have discy 
isophotes and often contain a rotating disc (e.g. 
Carter 1978, 1987; Davies et al. 1983; Bender et al. 



1988; Peletier et al. 1990; Jaffe et al. 1994). Core 
galaxies also tend to be more radio loud and have 
a greater soft X-ray flux (e.g. Ellis & O'Sullivan 
2006; Pellegrini 2010; Richings, Uttley & Kording 
2011, and references therein). 

It was, in part, from a diagram of central sur- 
face brightness versus magnitude that led Kor- 
mendy (1985b, his figure 3) to advocate a sepa- 
ration of dwarf and normal elliptical galaxies at 
Mb = —18 mag. However, as noted by Gra- 
ham & Guzman (2003), his sample was missing 
the bridging population near Mb = — 18 ± 1 mag. 
Excluding galaxies of this magnitude from Fig- 
ure [7^ would also result in two apparently dis- 
joint relations nearly at right angles to each other. 
It is therefore easy to understand how one may 
quickly reach the wrong conclusion from an incom- 
plete diagram. Although, Strom & Strom (1978, 
their figure 8; see also Binggeli et al. 1984) had 
already revealed that a linear relation exists be- 
tween magnitude and central surface brightness 
from —18.4 < My < —21.6 mag, spanning the 
magnitude gap in Kormendy (1985b). Nonethe- 
less, Faber & Lin (1983) had just observed that 
three of their six dwarf elliptical galaxies had near- 
exponential light profiles, leading them to specu- 
late that dEs are more closely related to "expo- 
nential systems" than (tidally truncated) ellipti- 
cal galaxies, and Wirth & Gallagher (1984, see 
also Michard 1979) had also just advocated a di- 
vision between exponential-dwarf and i? 1 / 4 -giant 
elliptical galaxies. 

To further confound matters, Kormendy (1985b) 
had the slope wrong for the distribution of dwarf 
elliptical galaxies in his M—fiQ diagram, which 
had two consequences. First, the bright end of his 
dwarf elliptical galaxy distribution did not point 
towards the faint-end of his luminous elliptical 
galaxy distribution, and thus there was no sug- 
gestion of a connection. This misrepresentation is 
unfortunately still propagated today (e.g. Tolstoy 
et al. 2009, their figure 1), although Kormendy et 
al. (2009) have now corrected this. Second, two 
points representing flattened disc galaxies were 
added at the bright end of the mis-aligned dwarf 
elliptical galaxy sequence by Kormendy (1985b), 
implying a connection between dwarf elliptical 
galaxies and disc galaxies rather than ordinary 
elliptical galaxies. 

A decade later, the Astronomy and Astro- 



21 



physics Review paper by Ferguson & Binggeli 
(1994; their Figure 3) had a big question mark 
as to "how" and indeed "if" dwarf and ordinary 
elliptical galaxies might connect in this diagram. 
When galaxies spanning the gap in Kormendy's 
(1985b) analysis were included by Faber et al. 
(1997, their figure 4c), and shown to follow a 
trend consistent with the relation from Strom & 
Strom (1978) and Binggeli et al. (1984), in which 
the central surface brightness became fainter with 
decreasing galaxy luminosity Faber et al. (1997) 
suggested that this behavior in their data was 
spurious and due to limited resolution — such 
was the belief in a discontinuity separating dwarf 
and ordinary elliptical galaxies. At the same time, 
Jerjen & Binggeli (1997) argued exactly the op- 
posite, suggesting that it was instead the "core" 
galaxies which had been displaced in the M-fi 
diagram from a linear M-fj,Q relation, rather than 
wrong central surface brightness measurements 
for the faint (non-dwarf) elliptical galaxies. As 
Graham & Guzman (2003) and Graham (2004) 
later explained, in terms of a ~0.1 percent central 
mass deficit relative to the outer Sersic profile in 
galaxies brighter than Mb ~ —20.5 mag, Jerjen & 
Binggeli (1997) were right, supporting the views 
expressed by Binggeli et al. (1984) on a continu- 
ity between dwarf elliptical and ordinary elliptical 
galaxies across the alleged divide at Mb ~ —18 
mag. 

3.3.2. Luminosity- colour relation 

Additional support for the dry merging scenario 
at the high-mass end is the flattening of the colour- 
magnitude relation above 0.5-1 x 10 11 M Q . While 
low luminosity, low Sersic index, elliptical galax- 
ies are bluer than bright elliptical galaxies (e.g. 
de Vaucouleurs 1961; Webb 1964; Sandage 1972; 
Caldwell & Bothun 1987), the brightest galaxies 
have the same colour as each other. This flatten- 
ing in the colour-magnitude relation was noted by 
Trcmonti et al. (2004) and is evident in Baldry et 
al. (2004, their Figure 9), Ferrarese et al. (2006a, 
their Figure 123), Boselli et al. (2008, their Fig- 
ure 7) and even Metcalfe, Godwin & Peach (1994). 
These observations help alleviate past tension with 
semi-analytic models that had predicted a rela- 
tively flat colour-magnitude relation for bright el- 
liptical galaxies (e.g. Cole et al. 2000). Previously, 
based on what was thought to be a linear colour- 



magnitude relation, Bernardi et al. (2007) had 
written that "if BCGs formed from dry mergers, 
then BCG progenitors must have been red for their 
magnitudes, suggesting that they hosted older 
stellar populations than is typical for their lumi- 
nosities". However, the flattening in the colour- 
magnitude relation has since been recognised in 
yet more data sets (e.g. Skelton, Bell & Somerville 
2009; Jimenez et al. 2011) although it should per- 
haps be noted that Skelton et al. reported the tran- 
sition at Mr = — 21 mag, i.e. ~1 mag fainter. 

In passing it is noted that the relation be- 
tween luminosity and supermassive black hole 
mass (Marconi & Hunt 2003; McLure & Dun- 
lop 2004) was found to be, after several refine- 
ments by Graham (2007), a linear one-to-one rela- 
tion for black hole masses predominantly greater 
than 10 8 M Q — consistent with the concept of dry 
galaxy merging at this high-mass end. 

3.3.3. Dynamics 

From a sample of 13 early-type galaxies, plus 
one spiral galaxy, Minkowski (1962) noted that a 
"correlation between velocity dispersion and [lu- 
minosity] exists, but it is poor". He wrote that 
"it seems important to extend the observations to 
more objects, especially at low and medium ab- 
solute magnitudes" . This was done by Morton & 
Chevalier (1973) who noted the same "continuous 
distribution of dispersions from 60 km/s for M32 
to 490 km/s for M87" but also did not attempt to 
quantify this trend. It was Faber & Jackson (1976) 
who, with improved data and a larger sample of 
25 galaxies, were the first to quantify Minkowski's 
relation and discovered that L cx cr 4 for their data 
set. This result has proved extremely popular and 
is known as the Faber- Jackson relation. Not long 
after this, Schechter (1980) and Malumuth & Kir- 
shncr (1981) revealed that the luminous elliptical 
galaxies followed a relation with an exponent of ^5 
rather than 4. At the same time, Tonry (1981) re- 
vealed that expanding the sample to include more 
faint elliptical galaxies results in an exponent of 
-3. This led Binney (1982) to write that "proba- 
bly the correlation cannot be adequately fitted by 
a single power law over the full range of absolute 
magnitudes" and Farouki et al. (1983) wrote that 
"the data suggests the presence of curvature in 
the L — a relation". Davies et al. (1983), and later 
Held et al. (1992), revealed that the dwarf ellipti- 
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cal galaxies followed a relation with an exponent 
of ~2, which explains why Tonry (1981) had found 
a slope of ~3 when including dwarf and ordinary 
elliptical galaxies. The relation found by Davies 
et al., with a slope of ~2, has recently been ob- 
served by de Rijcke et al. (2005) and the curved or 
possibly broken L — a distribution has been inter- 
preted by Matkovic & Guzman (2005) as a change 
in slope at —20.5 B-mag (see also Evstigneeva et 
al. 2007) in agreement with Davies et al. (1983). 

In spite of all the above work, there is a huge 
body of literature today which appears unaware 
that the L — a relation is curved or broken. Simu- 
lations of galaxies which succeed in producing the 
linear Faber- Jackson relation, L cx <r 4 , have actu- 
ally failed to produce the full distribution of dy- 
namics seen in real elliptical galaxies as a function 
of magnitude. 

4. Disc Galaxy Light Profiles 

Due to contractual agreements with Springer 
Publishing, the full version of this paper can not 
be posted to astro-ph. All of Section m plus the 
full version of section 13.3.11 will be available from 
the Springer Publication. 

4.1. The Bulge-Disc Decomposition 

4.2. Dust and inclination corrections 

4.3. Pseudobulges 

4-3.1. Sersic index 
4-3.2. Rotation 
4.3.3. Ages 
4-3-4- Scaling relations 

4.4. Bulgeless galaxies 

4.5. Barred galaxies 



5. Summary 

We have reviewed the progress over the last cen- 
tury in modelling the distribution of stars in ellip- 
tical galaxies, plus the bulges of lenticular and spi- 
ral galaxies and their surrounding discs. A num- 
ber of nearly forgotten or poorly recognised ref- 
erences have been identified. The universality, or 



at least versatility, of Sersic's R}l n model to de- 
scribe bulges (section 14. ip and elliptical galaxies 
(section 12.11) extends to the stellar halos of cD 
galaxies (section 12. 4[) and simulated dark matter 
halos fsection [2.1.ip . 

Dwarf and ordinary elliptical galaxies were 
shown in section 13.11 to be united by two con- 
tinuous linear relations between absolute mag- 
nitude and a) the stellar concentration quanti- 
fied through Sersic's (1963) R}l n shape parameter 
(section [3 . 1 ,2[) . and b) the central surface bright- 
ness, which is also related to the central density 
( section ^. l.lj) . As discussed in section [3~3l a break 
in the latter relation at Mb ~ —20.5 mag signals 
the onset of partially depleted cores relative to the 
outer Sersic profile in luminous elliptical galaxies. 
Additional scaling relations are also noted to show 
a change in character at this magnitude, which 
may denote the onset of dry galaxy merging. 

The identification of depleted galaxy cores and 
excess nuclear light relative to the outer Sersic 
profile was discussed in sections 12.21 12.2.11 and 
12.31 After accounting for these features, it was 
revealed how the above two linear relations result 
in curved scaling relations involving effective half 
light radii and effective surface brightness (sec- 
tion HO]), Specifically, the M-R c , M-fi e , M-(/x) e , 
/i e -R e , (ii} e —R e and n—R e relations are non-linear. 
These continuous curved relations exist because 
elliptical galaxies do not have a universal profile 
shape, such as an i? 1 / 4 profile, but instead a range 
of profile shapes that vary smoothly with absolute 
magnitude. Without an appreciation of the origin 
of these curved relations, they had in the past been 
heralded as evidence for a dichotomy between faint 
and bright elliptical galaxies. Numerical simula- 
tions and semi-numerical models which try to re- 
produce the full elliptical galaxy sequence must be 
able to reproduce these non-linear relations. This 
will likely require physical processes which work 
in tandem, albeit to different degrees over differ- 
ent mass ranges, to produce a continuum of galaxy 
properties that scale with mass while adhering to 
the linear M-n and M-^lq relations (subject to 
core-formation) . 

[snip] 

The upcoming 2.6 m VLT Survey Telescope 
(VST, Arnaboldi et al. 1998; Capaccioli et al. 
2005), plus the 4x1.8 m Pan-STARRS array 
(Kaiser et al. 2002), the 4 m Visible and Infrared 
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Survey Telescope for Astronomy (VISTA, Emer- 
son et al. 2004) and the 8.4 m Large Synoptic Sur- 
vey Telescope (LSST, Tyson 2001) are expected to 
deliver sub-arcsecond, deep and wide field-of-view 
imaging covering thousands of resolvable galax- 
ies. By pushing down the luminosity function 
into the dwarf galaxy regime, and through the 
application of improved galaxy parameterisation 
methods which allow for structural non-homology 
and the 2- or 3-component nature of disc galax- 
ies, both statistical and systematic errors will be 
reduced. This will undoubtedly provide improved 
constraints on galaxy scaling relations and, in 
turn, a fuller understanding of galaxy evolution. 

This research was supported under the Aus- 
tralian Research Council's Discovery Projects 
funding scheme (DP110103509). This research 
has made use of NASA's Astrophysics Data 
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NASA/IPAC Extragalactic Database (NED). 
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